Shear viscosity in late time of hydrodynamic evolution in AdS/CFT duality 
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We investigate the shear viscosity r\ and the entropy density s of strongly coupled TV = 4 su- 
per Yang-Mills (SYM) plasma in late time of hydrodynamic evolution with AdS/CFT duality and 
Bjorken scaling. We use correlation function method proposed by Kovtun, Son and Starinets. We 
obtain the metric in a proper time dependent AdSs space through holographic renormalization, 
whose boundary condition is given by energy-momentum tensor of the plasma in 2+1 dimension 
with transverse expansion or radial flow. With the metric we compute r\ and s of fluids in 1 + 1 and 
2+1 dimension without and with radial flow. We find the ratio 77/s = l/(4ff) in 1+1 dimension 
consistent with the Kovtun-Son-Starinets bound if next-to-leading terms in proper time are included 
in the equation of motion for metric perturbations. For 2+1 dimension the result is unchanged in 
the leading order of transverse rapidity. 

PACS numbers: 12.38Mh,11.25Tq 



I. INTRODUCTION 

Quantum chromodynamics (QCD) tells us that quarks and gluons are confined inside hadrons in vacuum or ground 
state. No free quarks and gluons are detected in any hadronic collisions. Finite temperature theory of QCD predicts 
that the deconfinement of quarks and gluons can be reached at high temperatures and densities and a new state of 
matter, the quark gluon plasma (QGP), can be created According to lattice QCD calculations the confinement - 
deconfinement phase transition takes place at temperature of about 170 MeV in three flavor case [2]. In early 1970s 
Lee and Greiner and his collaborators proposed that the deconfinement could be realized by smashing two heavy ions 
together at ultrarelativistic energies [jj. This is the only way of generating QGP in laboratories. The minimum energy 
density for deconfinement to take place is about an order of magnitude higher than normal nuclear matter density. The 
current energy frontier of heavy ion collisions is the Relativistic Heavy Ion Collider (RHIC) at Brookhaven National 
Laboratory (BNL), which has been running since the summer of 2000. A lot of evidences imply that the QGP at 
RHIC is near perfect fluid (strongly coupled QGP or sQGP), which is contrary to the conventional concept of QGP 
as a weakly interacting gas of quarks and gluons, see e.g. (author?) \4,m- 

While perturbation can be used to deal with weak coupling problems, there are a lot of difficulties in describing a 
strongly coupled system due to its non-linear feature and failure of perturbative methods. In recent years a promising 
method to deal with strong coupling problems in gauge theory is through string/gauge duality or precisely AdS/CFT 
duality, where AdS and CFT are abbreviations for anti-de Sitter space and conformal field theory respectively. The 
duality was first proposed by Maldacena and many others @, 0, Q, where a strongly coupled conformal gauge field 
theory corresponds to weekly coupled closed strings at 't Hooft limit. So properties of a weekly coupled string can 
provide information for strongly coupled gauge theory. 

Policastro, Son and Starinets first used AdS/CFT duality to calculate transport coefficients of strongly coupled 
dense matter @. Kovtun, and Son and Starinets derived a lower bound (called KSS bound) for the ratio of shear 
viscosity to entropy density, 77/s > 1/(47t), consistent to the RHIC data (author?) [rJ. These works open an avenue 
to study properties of stron gly coupled matter followed by a lot of developments in hydrodynamic properties of 
sQGP ( author ?) fill Ha Il4l|. Recently the jet quenching effect has been investigated with AdS/CFT techniques 
(author?) [13 [la. 17| . The sceening length of the heavy quark potential can also be described from the AdS/CFT 
duality (author?) HI, [H, H3, HH , which is closely related to the J/ip suppression in hot medium. 

The KSS bound is valid for any conformal fluids. Altough it is debatable that the sQGP is a conformal fluid, the 
ratio 77/s is fixed to be the KSS bound in all periods of the sQGP evolution in some hydrodynamic simulations [22L l23l| . 
However the sQGP formed in ultra-relativistic heavy ion collisions has particular geometric configurations such as 
boost invariance [13| etc.. Janik and Peschanski first considered a proper time dependent AdSz metric with boost 
invariance through holographic renormalization in late time of fluid evolution (author?) [25|, l26|, l27j . The method has 
been extended to sQGP with shear viscosity [2H |29|. where the shear viscosity appears in energy- momentum tensor 
in CFT and finally enter in the AdS$ metric. The early time behavior of sQGP has also been addressed following 
similar method and can then join the late time one to provide a global picture for time evolution in heavy ion collision 
(author?) [3(J [Sl|. The exact solution to the gravity dual of 1+1 dimensional sQGP with Bjorken scaling has been 
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found by Kajantie, Louko and Tahkokallio, which is helpful to understand what happens in 1+4 dimension in the real 
world (author?) [331 - The proper time dependent AdS$ metric without shear viscosity does not necessarily mean 
that the Af = 4 super Yang-Mills (SYM) plasma is non- viscous, noting that the shear viscosity with AdS/CFT duality 
is obtained by a D3 black AdS metric in Ref . [33l . l34T | . If we treat the proper time dependent AdS$ metric without 
viscosity as an extension of standard AdS metric we can build up a connection between the sQGP as an ideal fluid to 
the SYM plasma, then the ratio r]/s may not be zero. 

In this paper, we investigate the 2+1-dimensional expansion (also including 1+1-dimensional expansion) of the 
sQGP in late time and obtain the proper time dependent AdS 5 metric through holographic renormalization. The 
metric dual to the sQGP fluid can be obtained by solving the Einstein equation with the stress tensor of sQGP as 
boundary value of the metric if transverse expansion is small and can be treated as a perturbation. The metric has 
off-diagonal elements proportional to transverse rapidity. We compute with the metric the shear viscosity in 1+1 and 
2+1 dimension through Kubo formula with AdS/CFT duality. In 1+1 dimension with transverse expansion turned 
off, the shear viscosity is found to be vanishing if only the leading terms in proper time are considered in the Einstein 
equation or the equation of motion for perturbations to the metric. Only including the next-to-leading terms can 
one get the ratio r\js = l/(47r), indicating the shear viscosity is a higher order effect in late time behavior. In 2+1 
dimension the result is unchanged in the leading order of transverse rapidity. 

The structure of the paper is as follows. In Sec. [U we obtain the proper time dependent AdS$ metric corresponding 
to the sQGP in 2+1 dimension. In Sec. Hill we derive the effective action for perturbations in the AdS metric. We 
compute the ratio 77/s in Sec. IIV Al and IIV Bl for 1+1 and 2+1 dimensional cases respectively. We make summary 
and conclusions in Sec. El In the paper we take the signature (— , +, +, +, +) for the metric. 



II. 2+1 HYDRODYNAMICS WITH PROPER TIME DEPENDENT AdS s METRIC 



A. 2+1 hydrodynamics for ideal fluids 



In order to describe the hydrodynamic evolution of an ultrarelativistic system with transverse expansion, e.g. in 
heavy ion collisions, we use cylindrical coordinates, 



(r,y,r, 0), 



(1) 



where y is the rapidity, r the transverse proper time, r the radius in transverse plane, and 8 the azimuthal angle 



around longitudinal direction. Note that the full proper time 77 is related to r by rj 



r 2 . Hereafter the proper 



r 2 dy 2 + dr 2 



r 2 d0 2 . 



time means the transverse one if not explicitly stated. The line element is given by 

ds 2 = -dr 2 

which gives the metric gjS = diag(— l,r 2 , l,r 2 ). The velocity four-vector can be parametrized as 

dX>* 1 



(2) 



dr f ^l-R2( T ,r)/i 
[cosh a(r, r), 0, sinh a(r, r), 0] 



1,0,^,0 



(3) 



where R is a function of (r, r) and a is the radial rapidity in transverse direction. The radial velocity is defined by 
u r = sinha. One can verify that the velocity satisfies = — 1. 

For ideal fluid with the equation of state p = 3P where p and P are energy density and pressure respectively, the 
energy-momentum tensor is 

V - (p + P)u^ + Pgl°J(x) 



( — | + |pcosh 2 a — |p sinh a cosh a \ 
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(4) 



One can verify that satisfies conformal invariance = 0. The energy- momentum conservation equations read 



(5) 
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Figure 1: The radial velocity u r = sinha as functions of r and r. 
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Figure 2: The radial velocity u r = sinha as functions of r and r. 



where are Christoffel symbols for the metric in |f5|). The above equations lead to 



d T In p 



2 2 cosh 2 a 



t 2 cosh 



1 



d T lnR = l -*tL- 

T 



2 + 2(l-d r R) cosh 2 a 
2 cosh 2 a + 1 



(6) 



Given the initial condition i?(ro,r) = ^tqt where £ is set to 0.05 fm" 1 , Eq. ((6]) can be solved numerically whose 
results are shown in Fig. (Q]) ([2]) and are comparable with those in Ref. [35| . The results show that the radial velocity 
is linearly proportional to the radial distance, i.e. u r ~ £r, where £ is a small number. One can see in the figure 
that u r rises sharply in a very short early time and gradually falls with increasing r. If the expansion in transverse 
direction is negligible compared to that in longitudinal direction the energy density damps in the same way as in the 
case of 1+1 dimension, i.e. p ~ r~ 4 / 3 . 
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B. General setup of holographic model for 2+1 hydrodynamics 

In this section we briefly introduce the general idea of the holographic model for relativistic hydrodynamics. The 
metric in AdS space can in general be written as (36j, 

ds 2 = 9MN dX M dX N = \ [g^{x,z)dx il dx u + dz 2 ] , (7) 

where X M = (x^,z) and g,j, v {x,z) is the metric tensor in the 4-dimensional space ([I]) and generally depends on z. 
Normally one assumes the following form of 7j^ v {x, z), 

g^(x, z)dx {i dx v = -Adr 2 + Br 2 dy 2 + C (dr 2 + r 2 d0 2 ) + 2Ddrdr, (8) 

where we have introduced the off-diagonal term proportional to drdr to account for transverse expansion or radial 
flow implied by energy- momentum tensor in |(4]). Note that the off-diagonal term related to transverse expansion is 
absent in the metric (J2j) - One can write the coefficients A, B, C and D in exponential forms, 

A = e w , B = e\ C = e c , D = e d . (9) 

In general these coefficients are functions of x and z. A standard way of holographic renormalization [37I . l38| is to 
expand the metric g^(x, z) or these coefficients in the fifth coordinate z, 

g llv {x,z) = Y,^ n 9 i ^ ) {x), (10) 

71=0 

where gjw(x) is just the metric (J2J) in 4-dimensional space. The second order term gj?J(x) can be proved to be 
vanishing. The fourth order term gffi (x) is given by the energy-momentum tensor, 

sit'w^w, (ii) 

where we set factor ^ = 1 (N is the number of colors) and we will finally restore it. Given by these boundary 
conditions one can solve the AdS metric J7]) from the Einstein equation with the cosmological constant A = —6 which 
we will discuss about in Sec. llIIII) . 

Rmn — -^qmnR + QgMN — 0, (12) 

where Rmn and R are the curvature tensor and scalar in AdS space with the metric gMN given in (0. 



C. Holographic model in late time 

The introduction of the r dependence of the metric (|22|) makes solving the Einstein equation a formidable task even 
with the method of Janik and Peschanski (23, l28l. |29| . So the simplification is necessary. In this section we consider 
a simplified problem with perturbation in transverse direction based on the Janik-Peschanski method. 

Supposing a is small in Eq. ((4]) we can make an expansion of T M „ in a, 
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(13) 



In this expansion, we see that off-diagonal elements are proportional to odd power of a, while modifications of diagonal 
elements are proportional to even power of a. In leading order the transverse velocity vanishes and we reproduce the 
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results in Ref. [2a.l2a|. In next-to-leading order 0(a) the transverse velocity does not vanish while the energy density 
and pressure do not change, whose modifications are in 0(a 2 ) or higher. 

If we assume that there is no linear terms of a in diagonal components of the AdS metric §m n, we get the following 
result by solving the Einstein equation (fl2|) . 



9mn — —x 

z A 



( 

1+a 

T 2 (l + fl) 

4„,(l-a) 2 n 



4 (1-a) 
3 a 1+a 



V 



o o \ 



1 + a 

r 2 (l + a) 

1 / 



(14) 



with a = ^ 4/3 and po is a constant energy density. 

In the case without off-diagonal components in the metric (fll|. i.e. a = 0, the line element in central rapidity 
region (not for non-central rapidity region) has the form, 



, , (1 - ZVZ% Y dT 2 1 + Z A /Z% , 2 , n , 2 2 dz 2 

where the black hole locates at 

*, = (?)-">*, 

which depends on r. See Appendix (jAj) for the derivation of the metric (fl5l) . Changing variables, 



ZH ^ Z H = 



1 



V2' 



(15) 



(16) 



(17) 



we can verify that the line element ifTB]) is the standard one in D3 black AdS$ space if z is constant, see Sec. 



ds 2 



dT 2 



dx 2 



dz 



r/2 



1 - Z 4 /^f Z 2 ' 



(18) 



where dx 2 = T 2 dy 2 + dr 2 + r 2 ^ 2 and zh is the location of the black hole horizon, the cosmology constant for AdSd+i 
is 1381 



A = - 



d(d-l) 



G. 



and the Hawking temperature is determined by the behavior of the metric near the horizon, 

1 V2 



T, 



H 



itzh kzh 



(19) 



(20) 



Hereafter we call the proper time dependent AdS§ metric the AdS* one in (fl4"ll with transverse expansion or in 
IH) without transverse expansion. 



III. EFFECTIVE ACTION FOR PERTURBATIONS IN AdS METRIC 



The low energy limit of type IIB string theory in AdS§ x S 5 space can be approximated by a five-dimensional theory 
(see, e.g. (author?) j34[) whose action is 



iv 2 r , 

ho^^^J d 5 x(U 5D -2A + ---), 



(21) 



where TZ^d is the scalar curvature in AdS§, A the cosmological constant, and R the radius of S 5 . The line element in 
AdS§ space without black hole is 



r 2 R 2 
ds 2 = -=^{-dt 2 + dx 2 ) + -^dr 2 . 
R z r z 



(22) 
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Using the variable z = R 2 /r, the above line element can be written as 

I 2 

7 2 

The line element with a black hole has the form 



p2 

2 / i 2 2 2 \ 

rfs = — =-(— at + dx 4- ' 

z z 



(23) 



d S 2 = ^ 



{1- z 4 /z 4 H )dt 2 +dx 2 + 
^ K [-(l-^ 2 )dt 2 + rfx 2 ] + 



1 



1-^/4 
i? 2 



4m 2 (1 - m 2 ) 



cfe 2 
du 2 , 



where it = z 2 /z 2 1 and = l/(wT). 

The AdS/CFT correspondence can be expressed by 

Z iD [J] = e iS ^ 



J(x)=<j>(x,z=Q) 



(24) 



(25) 



where Zhe>[J] is the partition function of CFT in 4-dimensional Minkowski space and Ssd[4>] the action in AdS$ given 
by (pi]) for the classical bulk field <j)(x, z), J(x) is the source in 4-dimension taking the value of the bulk field on the 
boundary z = in AdS§. The partition function Z±e>[J] is a functional of J{x), 



Z±d[J] = J [rf^exp I^cftM d A xJO 



(26) 



where <p denotes the CFT fields, 0{<p] is the operator coupling to J(x), and Scft[<p] is the CFT action. Following the 
AdS/CFT correspondence the Green functions of operators O in CFT can be evaluated in terms of S^z?^] in AdSs, 
for example, a two-point Green function is 



{TO{ Xl )0(x 2 )) AD = 



S 2 S 5D [ 



/4D iStfaWfa) 
We note that the stress tensor is coupled with the metric, 



J(x)=cj>{x,z=Q) 



T„ 



2 SScft 



(27) 



(28) 



fJ 



where can be regarded as an operator in CFT and the metric g^ v in 4-dimension as its source. In order to obtain 
the shear viscosity 77 through the Kubo formula, 



n 



lim -ImGf 212 (w,0) 



we need to know the Green function of Ti 2 at low energy limit, 



G 



12,12 



(w, k) = -i 



d A xe lkx 



-e(i)<[T 12 (z),T 12 (0)]>. 



(29) 



(30) 



To obtain Gf 2 12 ( w jk) from the gravitational dual we should know the action S^olg^u] in AdS§ as a functional of 
g^u with the boundary value g^ v {z = 0) — > g^. For this purpose we consider a variation in the AdS§ metric, 



9nv — g$v + h^ v . We can choose for simplicity h z ^ = h^ z = 
element in 11241) then becomes 

_2t^2 p2 

ds 2 = [-(1 - u 2 )dt 2 + dx 2 ] 



and that h^ v depend on t, x 3 and z (or u). The line 
R 2 



4u 2 (l - u 2 ) 



du + haudx^dx 1 ' . 



(31) 



We can expand the action lj2T|) in terms of h^u, 

iV 2 
8^ 

8^ 



I5D 



d 5 x^(U^ - -g^K + Ag^W 

d 5 xj— g {nf u - l -g^n w + \ g $)hr 



m2 p 1 1 



8tt 2 

r2 



;y2 p 1 

^ I d 5 x^9(Sn^--g^5n + Ah^)h^, 



(32) 
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where we have used the fact that the first term of the second equality is vanishing from the Einstein equation for the 
unperturbed metric g'f}. We have also used TZq = —20 and A = —6. Here we have chosen the S 5 radius R — 1. We 
consider the components h\2 — h 2 i which can be verified to decouple from others. For the metric lf3lj) . the action 
<f32|) can be simplified as 

hD ~8^J d " x ^9^ R ^ + 4h 12 )h 12 = J d 5 xV=9 i-^d^d^ + ■■■), (33) 

with (f) — h\ = g xl h\2 and h 12 = —g 11 g 22 h\2- The retarded two-point Green functions of the CFT stress tensor can 
then be derived from the quadratic terms in I$n as above by using Eq. (|27|) , The shear viscosity is then obtained 
through the Kubo formula (|29)l . 



IV. SHEAR VISCOSITY AND ENTROPY DENSITY IN AdS* 

In this section we will calculate the shear viscosity in AdSf following the procedure of AdS / CFT duality presented 
in the previous section. Also we will compute the entropy density in AdSf . We consider fluid evolutions in 1+1/2+1 
dimension without/with transverse expansion or radial flow. 



A. 1 + 1 dimension without transverse expansion 

Following Janik and Peschanski's solution [HI in late time, the line element in AdSf in 1+1 dimension can be 
written as 

ds 2 = i \-^-^dr 2 + (1 + a)[r 2 dy 2 + dx\ + dx 2 ]\ + (34) 
z l \ 1 + a ) z z 

4 

where a = ^ 4/3 . Note that the derivative of a is non-zero, i.e. d T . z a ^ in computing the Ricci tensor and the 
scalar curvature. Finally we will keep a constant a while taking the limit r — > oo when doing power counting in 
t. In adopting the above measure we can verify after a lenghy but straightforward algebra that the Ricci tensor 
IZ^ = — Ag^u and IZ^d = —20 upto higher order terms in negative powers of r. The determinant of the metric is 

= A-(l — a 2 )r. Taking a perturbation in the metric <? M „ = g$ + h^ v with hyi — hnij, y, z) and h^v = for all 



other indices fx, v, we obtain the quadratic terms of the effective action up to higher order contributions in negative 
powers of r, 

ho « -t^t / tfxyFlgrd^h (35) 



16tt 2 j 

where <p = h\ = g ll h\2, d 5 x = drdydxidx 2 dz and /i = T,y, 1, 2, z. The equation of motion for 4> is then 

d^V^gg^d^) = 0. (36) 

We assume a factorized form for <j>, 

, 2 J exp(-iwrcoshy + ip 3 r sinhy)^o(w,p 3 )/ P (^) (37) 

where p 3 is the momentum along the third axis (we do not distinguish subscript or superscript for this mometum). 
Inserting the above into Eq. f36|) and taking the limit r — > 00 while keeping v or a constant, the leading order 
contribution of 0(t°) in Eq. |(4|) gives 

° = (1 ~ ° 2) ^ + fp ^1 -"a ^ C ° Sh2 V ~ " Q )(- LJsinh ^ +P3 coshy) 2 . (38) 
We consider the central rapidity y = and static limit p 3 — 0, the above equation becomes 



d 2 f P , (l+o) 



+ ;r 1 ,> 2 / p = 0. (39) 



dz 2 (I -a) 
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The reason for choosing central rapidity is that the metric ([34)1 can only be treated as an extension of the standard 
AdS one at central rapidity, see Appendix [A] Considering the solution near the horizon at aa = 1 and changing 
variables, 



3 
A 

2V2 V 3 



^)~ 17 V/V (40) 



Eq. (|39|) can be rewritten 



The solution to Eq. (|4Tj) can be found, 



/p(z) = (1 - z ')l/2(l±Vl-4-' 2 ) M (1 _ z')l/2±W (42) 

where we have approximated — 4o/ 2 « i2o/ at the limit r — > oo, since cj' ~ r 1 / 3 3> 1. Using the incoming wave 
solution corresponding to the positive sign in (|42| and substituting it back into the action (|35| , the boundary term at 
z = zh or the black hole horizon at an — 1, the retarded Green function 12 (w) at the static limt can be obtained, 
which is vanishing due to the presence of the real part \ in the exponent of the solution lj42j) . This can be seen from 
the fact that there is a factor (1 — z') which is zero at the boundary, see Eq. IjBip in Appendix ([B|. Therefore we 
have shown that the shear viscosity is absent from the scaling solution in 1+1 dimension in the leading order. 
Picking up the next-to-leading order contribution of C^I/t 1 / 3 ) in Eq. l(36|) . Eq. I|39p becomes 

<Pf p 3 + 5a 2 df p (1 + a) 2 

d^""z(i- a 2) -di + jr^ u fp -°- (43) 

Changing variables as in Eq. l(40|) . the above equation near the horizon can be rewritten in the form 

d 2 f P I df p 1 , 2 



whose solutions are 



dz* 1 - z> dzJ ' (1 - Jf* ** = °> (44) 

/ ±p (z) = (l-z') W - (45) 
As derived in Eq. (|B2|I in Appendix [Bj the retarded Green function at static limit is obtained, 

Gf 2 , 12 K0) = -2F( W ) = -i ^^ ^/V 1 ^. (46) 
The shear viscosity per unit transverse area can be obtained from the Kubo formula (|29)l . 

, = ^i (47) 

where we have recovered the factor p — ► po (i^) 

The Hawking temperature is T = from Eq. (|20|) where Zjj is the horizon of the black hole given in Eq. (fTB")) . 
Then the initial energy density is a constant due to T ~ t -1 / 3 , 



7T \ 3A^ 2 

The entropy per unit rapidity and unit transverse area is given by jH, E 
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The entropy density is obtained, 

8 = 7 = \ N >" T ^ ( 5 °) 

where we see that the entropy density has an asymptotic behavior s ~ r _1 . From Eq. f47j) and (|50|l . we get the 
well-known value, 

- = f- ( 51 ) 

S 47T V ^ 

B. 2+1 dimension with transverse expansion 



The metric in 2+1 dimension with radial flow is given in l|14p . in order to calculate the shear viscosity we have 
to introduce a perturbation to the background metric. It is then convenient to explicitly use rectangular transverse 
coordinates (2:1,3:2) instead of cylindrical ones (r, 6). The corresponding line element is 



If (1-a) 2 ,„ w , , , , . 8 (l-a) 2 l , , dz 2 

— < at + (1 + a){dx^ax^ + dxidxi) + —a — Xidxidt H — 

z z I 1 + a 3 1 + a r z 



where r — ^fxixl and the summation over i = 1, 2 is implied. We keep in mind that a is small. The metric determinant 
has the same form as in the 1+1 dimensional case in the leading order, 

^—g=± E (l-a 2 ) + 0(a 2 ). (53) 

Now we consider the perturbation hi2(t, X3, z) to the background metric (|52|) . Denoting 4> = h\(t, X3, z), we find the 
equation of motion for <fi, 

= ST^ + 4^= -— L = a /i (V=s^"ft,0). (54) 
9 

It can be verified that the entangled term ~ acf> appear in 0(t -1 / 3 ), while leading terms are of 0(i 2 / 3 ). So the 
transverse part is decoupled from (f> in the leading order. The shear viscosity and entropy density per unit transverse 
area are the same as in the 1+1 dimensional case, Eqs. l(47|) and ([50|l . 

We can also consider the perturbation along X2 axis (or equivalently X\ axis), i.e. 4> = h\(t, z, X2). To the leading 
order 0(r 2 / 3 ) we find the equation of motion for <f>, 

= 5K\ + 40 = --J—d^^dg^d^) + t^loc-^—^-dt^ (55) 



2V=9 3 1 + ar 

whose explicit form reads, 



We see that the transverse part enters the equation of motion in the leading order. We assume that the solution has 
the factorization form (f>(t,X2,z) = 4>a{p)fp{z)e~ tuJt+lp2X2 , then we derive from Eq. 1(56]) the differential equation for 

d 2 f p 3 + 5a df p 1 + a 2 i 2a 

tiy^)lTz + (i - ^ fp TT^3r p = °' (57) 

where we can expand the equation near the horzion z — * zu in Eq. (fl6|) . 



(58) 
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The solution is a linear combination of Bessel functions, 

f P i = J iu < i 3/2 (l-z')z H 
f P 2 = Y iw , i 3 / 2 (l-z')z H 



Since a is small, we can make expansion in a for the solutions f p i and f P 2. One can verify that the solution is a linear 
combination of f p and 




f-P = fp( iuJ ' -> 



r(i + «v) (1 + ia/)r\i + 



aw 
37 



(59) 



One sees that /_ p = /*. Following Eq. (|B3[) in Appendix [B] and steps in previous section, we get the same value as 
in the case of 1+1 dimension, i.e. rj/s — l/(47r). 



SUMMARY AND DISCUSSIONS 



We derive a time dependent metric dual to sQGP fluid in 2+1 dimension with radial flow in late time by holographic 
renormalization. It is difficult to obtain the exact solution to the Einstein equation with this metric, especially when 
the metric has off-diagonal components for radial flows. If transverse expansion is small and can be treated as a 
perturbation, the late time asymptotic solution, the metric in (fT4|) with off-diagonal elements, can be found by using 
v = z 4 /t 4 / 3 as a scaling and expansion parameter. With this metric we calculate the ratio rj/s of shear viscosity 77 to 
entropy density s for sQGP in SYM field theory with the KSS method. As a first attempt we consider 1+1 dimension 
with only longitudinal flow whose metric is diagonal. If we include only the leading order terms in the equation of 
motion for perturbations to the metric, the shear viscosity is vanishing, consistent to the assumption that the fluid is 
ideal. We reproduce KSS bound 1 / (47r) for r\j s if we pick up the next-to- leading order term in the equation of motion, 
indicating that the shear viscosity is a higher order effect. Our derivation is based on the Janik and Peschanski's 
method and is valid in late time, r — > 00. For intermidiate stage of hydrodynamic evolution the ratio is not necessarily 
l/(47r) for an expanding fluid, so our result is not trivial or obvious. We further show that the ratio for fluids in 2+1 
dimension in late time with transverse flow is the same as in 1+1 dimension in the leading order of transverse rapidity. 
We remember that the mean free path is Z m f p ~ l/(nav), where n ~ T 3 is the particle number density, a ~ g 2 T~ 2 (g 
is the coupling constant) the typical scattering cross section and v ~ 1 the typical velocity. The shear viscosity can 
then be estimated as 77 ~ pl m f P ~ g~ 2 T 3 , where p is the energy density. In comparsion with our result in Eq. (|iT|l . we 
have T) ~ r _1 , T ~ r -1 / 3 and g does not change with the proper time. This implies that if the coupling is as strong 
at the beginning as in late time of fluid expansion or hydrodynamic evolution does not influence the strength of the 
interaction. 

In 1+1 dimension one can introduce the shear viscosity of the next-to-leading order in the stress tensor in late time 
solution (author?) [l^, [H, ■ So an interesting attempt is to calculate the shear viscosity with the metric dual to 
the stress tensor with shear terms. We found that the correction to the shear viscosity is also of the next-to-leading 
order, rj/s = 1/(4tt) + O^" 1 " 3 ). 

The same analysis can also be applied to hydrodynamic behaviors of sQGP in early time, which is important to 
understand the initial state of QGP (author?) [3(3] ■ However early time behaviors of fluids show anisotropic feature 
and therefore more complicated. We will reserve it for a future investigation. 
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Appendix A: DERIVATION OF THE METRIC IN ([HD 

The D3 black AdS metric can be writen as 



cfac 2 dz 2 

2~' 



Z^ 



Z^ 



Changing variables to (t 7 x 3 ) 
ds 2 



where A = z / f°J and I? = (1 + z^/zq). The metric (IA2I) has the symmetry under y — > —y. In comparsion with 

Janik and Peschanski's time dependence metric, there is an off-diagonal part ~ drdy in the metric l|A2|) . So the metric 
(fT5| is only valid in central rapidity region around y ~ 0. 



l + zV4 

(r, y) through i = rcoshy, x 3 = rsinhy, the metric becomes 

(-Acosh 2 y + Bsinh 2 y)dr 2 + r 2 (-Asinh 2 y + Bcosh 2 y)dy 2 
+2t(—A + B) cosh y sinh ydrdy + ... 



(Al) 



(A2) 



Appendix B: EVALUATION OF RETARDED GREEN FUNCTION FROM SOLUTION (|42|) 



Substituting the solutions f± in (|42|) and |45|) back into the action IpTj) and keep the boundary term at z = Zh, we 
find 



N 2 f 4 d 

167T 2 / ^ XdZ ^[ v/ ~5 5ZZ ^( T ' y ' 2: ) 5;5 ^(' r ' y ' Z )] 



TV 2 
16V 2 

TV 2 
16tt 2 



1 - a 2 

d xt ^—(j>{T,y,z)d z (j)(r,y,z) 



drdydx 1 dx 2 r 



dwdp^duj' dp' 3 



N 2 



\ Q n 2 J d(TCOshy)d(Tsmhy)dx 1 dx 2 J — ^ — s - 

x exp[-i(w + u')TCOshy + i(p 3 + p 3 )r sinh y](f> (Lj\p 3 )(j) Q (uj,p 3 ) 



exp[— i(u> + uj')t cosh y + i(p 3 + p' 3 )t sinhy] 



dujdp 3 duj'dp 3 



1-q 2 

N 2 
16?r 2 

(2*)' 



-/ P <(2) 



dz 



,1,1 du>dp 3 . 1 d 2 , d/ p (; 



■(/>o(-uj, -p)<j) (uj,p)F(u). 



In the fourth equality we have changed the integral variables, drdyr — d(r cosh y)d(T sinh y). In the last we have 
assumed that 4>q{lo,p) = (2tt) 2 8{pi)5{p2)4>v{—u 1 —pz) and set the transverse area L 2 = 1. Now we evaluate F(u). For 
the solutions (1421. we have 



F(u) oc hm y 3 ^ z 



z^z H Z u 



4 z ^ 



lim (l-z')(l-z') 1/2 ^'f (!-*') 



- d 



dz 



-^ t- 4/3 (1/2 + z^) lim (1-z') 



(Bl) 
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while for the solutions H45D. we obtain 



N 2 1-a- 



For the solution (|59| . we get 



where we have used 

dfpjz) 
dz' 



16tt 2 

N 2 



-f- P {*) dJM 



dz 



T 2 U 3 W 



4 / 3 lim (1 - z')(l - z')- iu3 \l - zj J ~ x 

Z^ZH 

//'o\ 3 / 4 _ a 

8 ^ lyJ r • 



1()77- J // 

. V^iV 2 /£>n\ 3 / 4 



iV 2 l-a : 



16?r 2 

iV 2 8 
16?r 2 3 

V2N 2 



-f- P {*) dJM 



dz 



p T-^ 3 lim (1 - z')f- p {z) 



dfpjz) 
dz' 



Po\ 3 / 4 _i 



V 3r 



(B2) 



(B3) 



aw 



(_1)^74 2 - JW ' i( _ 3^^3^742-2-^' 2 

r(i + zw') + (i + ^')r(i + ~ z ^ Zff_ 37 



-2(1 - z')*" +2 



■2—iuj 



3r / (1 + iw')r(l + iw') 



[1] T. D. Lee and G. C. Wick, Phys. Rev. D 9, 2291 (1974). 

[2] F. Karsch, E. Laermann and A. Peikert, Phys. Lett. B 478, 447 (2000) [arXiv:hep-lat /0002003| . 

[3] J. Hofmann, H. Stocker, W. Scheid and W. Greiner, Bear Mountain Workshop, New York, Dec 1974. 

[4] M. Gyulassy, L. McLerran, Nucl. Phys. A750, 30(2005)^ 

[5] E. V. Shuryak, Nucl. Phys. A 750, 64 (2005) |arXiv:hep-ph/0405 066 1. 

[6] J. M. Maldacena, Adv. Theor. Math. Phys. 2, 231 (1998) [Int. J. Theor. Phys. 38, 1113 (1999 )] |arXiv:hep-th/9711200] . 
[7] S. S. Gubser, I. R. Klebanov and A. M. Polyakov, Phys. Lett. B 428, 105 (1998) [arXiv:hep-th/9 802109l. 
[8] E. Witten, Adv. Theor. Math. Phys. 2, 505 (1998) JarXiv:hep-th/9803131] . 

[9] G. Policastro, D. T. Son and A. O. Starinets, Phys. Rev. Lett. 87, 081601 (2001) |arX iv:hep-th/ 0104066|. 
[10] P. Kovtun, D. T. Son and A. O. Starinets, Phys. Rev. Lett. 9 4, 11 1601 (2005) |a rXiv:hep-th/0405231| . 
[11] A. Buchel and J. T. Liu, Phys. Rev. Lett. 93, 090602 (2004) [arXiv :hep-th /03 1 1 1 75 1 . ~ 
[12] K. Maeda, M. Natsuume and T. Okamura, Phys. Rev. D 73, 0660 13 (2006) [arXiv:hep-th/06020"To] . 
[13] M. Natsuume and T. Okamura, Phys. Rev. D 77, 066014 (2008) |arXiv:0712'2916l [hep-th"]]~ 

[14] R. Baier, P. Romatschke, D. T. Son, A. O. Starinets and M. A. Stephanov, JHEP 0804, 100 (2008) |arXiv:0712.245TI 
[hep-th]]. 

[15] H. Liu, K. Rajagopal and U. A. Wiedemann, Phys. Rev. Lett. 97, 182301 (2006) | arXiv:hep-ph/ 0605178l. 

[16] C. P. Herzog, A. Karch, P. Kovtun, C. Kozcaz and L. G. Yaffe, JHE P 0607, 013 (2006) arXiv:hep-th/060 5158| . 
[17] S. S. Gubser, Phys. Rev. D 74, 126005 (2006) [arXiv: hepTh/0605182] . 

[18] H. Liu, K. Rajagopal and U. A. Wiedemann, Phys. Rev. Lett. 98, 182301 (2007) |arXiv:hep-ph/0607062|. 

[19] K. Peeters, J. Sonnenschein and M. Zamaklar, Phys. Rev. D 74, 106008 (2006) |arXiv:hep- th/0606195P 

[20] D. Hou and H. c. Ren, JHEP 0801, 029 (2008) [arXiv:0710.2639 [hep-ph]]. 

[21] M. Li, Y. Zhou and P. Pu, JHEP 0810, 010 (2008) [arXiv:0805.1611 [hep-th]]. 

[22] H. Song and U. W. Hein z, larXlv:0712.3715l [nucl-th], 

[23] A. K. Chaudhuri. larXiv:0801.3"T80l [nucl-th], 

[24] J. D. Bjorken, Phys. Rev. D 27, 140 (1983). 

[25] R. A. Janik and R. Peschanski, Phys. Rev. D 73, 045013 (2006) |ar Xiv:hep-th/0512162| . 
[26] R. A. Janik, Phys. Rev. Lett. 98, 022302 (2007) [arXi v:hep-th/06 10 144| . 

[27] P. Benincasa, A. Buchel, M. P. Heller and R. A. Jan ik, Phys. Rev. D 77, 04 6006 (2008) |arXiv:0712. 20251 [hep-th]]. 

[28] S. Nakamura and S. J. Sin, JHEP 0609, 020 (2006) |arXiv:hep-th/0607123| 

[29] S. J. Sin, S. Nakamura and S. P. Kim, JHEP 0612, 075 (2006) |arXi v:hep-th/0610113| . 

[30] Y. V. Kovchegov and A. Taliotis, Phys. Rev. C 76, 014905 (2007) |arXiv:0705.1234f Thep-ph]]. 



13 



[31] J. L. Albacete, Y. V. Kovchegov and A. Taliotis, JHEP 0807, 100 (2008) [arXiv:0805.2927 [hep-th]]. 
[32] K. Kajantie, J. Louko and T. Tahkokallio, Phys. Rev. D 76, 106006 (2007) [arXiv:0705.1791 [hep-th]]. 
[33] D. T. Son and A. O. Starinets, JHEP 0209, 042 (2002) |arXiv:hep-th/020505l|. 
[34] D. T. Son and A. O. Starine ts, Ann. Rev. Nucl. Part . Sci. 57. 95 (2007) |arXiv:0704.0240l [hep-th]]. 
[35] P. F. Kolb and U. W. Heinz, [a7xiv:nucl-th/0305084| 

[36] D. Z. Fefferman and C. R. Graham, Conformal invariants, Elie Cartan et les Mathematiques d'aujourd'hui (Asterisque, 
1985). 

[37] A. Karch, A. O'Bannon and K. Skenderis, JHEP 0604, 015 (2006) arX iv:hep-th/0512125| . 
[38] K. Skenderis, Class. Quant. Grav. 19, 5849 (2002) [arXiv:hep-th/0209067| . ~ 



